clear;clc;
for b=1:995
    a=sqrt((b+7)^2-b^2);
    if(a==floor(a))
        fprintf('a= %i,b= %i,c= %i\n',a,b,b+7)
    end
end

fprintf('\n')
syms u b
a=7*u;
b=solve((a^2+b^2-(b+7)^2),b)
c=sqrt(a^2+b^2)

b=(5*u^2)/2 - 5/2;
c=(25*u^2 + ((5*u^2)/2 - 5/2)^2)^(1/2);

% for u=1:50
%     a=5*u;
%     b=(5*u^2)/2 - 5/2;
%     c=(25*u^2 + ((5*u^2)/2 - 5/2)^2)^(1/2);
%     if(c<=996)
%         fprintf('a= %i,b= %i,c= %i,   i= %i\n',a,b,c,u)
%     end
% end